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SUMMARY 

This  report  describes  applications  of  the  method  of 
characteristics  to  supersonic  flow  about  axially  symmetric, 
al,  open-nosed  bodies,  such  as  diffusers.  A  general  proce¬ 
dure  is  given  for  the  computation  of  velocity  distributions 
on  the  external  surface  of  axially  symmetrical  diffusers  of 
arbitrary  shape.  The  case  of  straight-line  diffusers  is  de¬ 
veloped  in  more  detail,  including  the  case  with  incidence. 
A  special  function  is  introduced  which  allows  a  quick  ana¬ 
lytical  solution  of  the  problem  of  straight-line  diffusers. 


Ill 
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1,  INTRODUCTION 


The  main  purpose  of  this  paper  is  the  discussion  of 
the  axially  symmetrical  supersonic  flow  past  open  nose  dlf 
fusers  (with  the  exception  of  the  last  section,  where  also 
the  solution  for  flow  with  an  incidence  angle  is  outlined). 

The  diffusers  are  assumed  to  have  a  sharp  lip  (see 
figure  1)  and  such  shape  that  the  flow  is  isentropic  and 
may  satisfy  .pertain  simplifying  assumptions  discussed  in 
section  2. 

In  discussing  the  flow  past  diffusers  only  the  exter~ 
nal  one  is  considered.  For  this  to  be  justified,  perfect 
matching  must  be  assumed;  i.e.  there  is  no  "spillage'* . 

The  x,y  reference  system  is  chosen  in  such  way 
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Fig.  1 

that  the  x  axis  coincides  with  the  symmetry  axis,  and  has 
the  same  sense  as  the  undisturbed  flow  velocity. 

The  local  surface  angles,  formed  by  the  tangent  of 
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the  generator  of  the  diffuser  outer  surface,  with  the  sym¬ 
metry  axis,  are  called  6  . 

All  the  conditions  at  the  nose  point  are  designated 
with  the  subscript  "o"j 

e.g.  y6 ,  0„  ,  u,  ,  etc. 

Notations  and  definitions  follow  those  normally  used. 
Part icularly 

W  -  represents  the  value  of  flow  velocity  at  any 
point  of  the  flow  field. 

Wo,  -  the  value  of  velocity  in  the  undisturbed  region. 

U,V  -  represent  the  horizontal  and  vertical  compo¬ 

nents  of  the  velocity  W,  respectively. 

u,v  -  represent  the  horizontal  and  vertical  compo¬ 
nents  of  the  perturbation  velocity. 

Due  to  coincidence  of  the  direction  of  W  <*,  and  the  x 
axis,  it  is  possible  to  form  the  following  relations: 

n„  =  wTO 
=  0 

U  =  Uw  +  u 

V  =  v 

The  symbol  |3  denotes  cot  p.  ,  where  jjl  is  the  Mach 
angle  j  M-denotes  the  Mach  number  (  fi  —  V M1- d  ’  )  .By  "Mach 
lines"  we  shall  understand  the  straight  lines,  tangent  to 
the  characteristic  lines.  Both  lines  (Mach  and  character¬ 
istic)  are,  of  course,  identical  in  the  undisturbed  flow 
region. 
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In  later  sections  the  "straight  line  -cliff ■user"  will 
be  separately  discussed.  By  this  name  'we  shall  '•understand 
those  diffusers  whose  outer  surface  forms  a  straight line 
in  an  axial  section  (see  figure  7)»  '•-■■■■ 

Here  is  a  short  outline  of  the- paper:'  - 

In  sections  2  to  8  the  simplified  characteristic  equa 
tions  according  to  the  Sauer -Heinz'-  method  are  adapted-  to 
diffusers.  ’  •  •"  •  1  ■  1  ' 

Nets  of  Mach  lines  and  equations  for  u and  V  are  es¬ 
tablished  for  each  point  of  net  which  allow  the  naldula-' 
tion  of  these  velocities  from  known  data  at  former  points 
of  net.  Next,  by  means  of  algebrical  transformations,  the 
v  -  components  of  pert  ion  velocity  are  eliminated, and, 
in  special  cases,  solutions  are  found  for  u  directly  from 
the  geometrical  position  of  the  points  ("the  first  Mach 
line",  see  section  5)» 

After  adapting  the  above  described  method  to  straight 
line  diffusers  (section  9),  the  possibility  is  outlined  of 
solving  this  problem  (i.e.  of  straight  line  diffusers)  by 
analytical  means;  this  is  facilitated  by  assuming  a  func¬ 
tion  cr  (x,  ©,|3)  ,  defined  by  equation  (16),  section  10, 

and  its  limit  <r  (x,[i^  for  9^*0  (section  11). 

The  next  sections  are  dedicated  to  finding  the  numer¬ 
ical  values  of  function  CT(x,  (i)  (section  12)  and  its 
analytical  form,  first  on  the  basis  of  formerly  described 
methods  (section  13),  and  next,  utilizing  operational  cal- 
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cuius  (sections  li|-l6).  There  are  obtained  an  exact  solu¬ 
tion-equation  (39)  and  approximate  solutions-equations(22), 
(kO)  and  (ill). 

In  sections  17-18  the  approximate  solution  of  the  u- 
velocity  distribution  is  given,  not'  only  for  the  diffuser 
surface,  but  for  the  entire  space  around  it. 

Finally,  in  the  last  section,  the  solution  of  flow 
with  an  incidence  angle  is  outlined. 

Graphs,  ilustrating  all  more  important  results,  are 
enclosed. 
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Part  One  -  GENERAL  METHODS 

Z.  THE  DERIVATION  OF  THE  SAUER-HEINZ  EQUATION 

The  considerations  of  this  section  will  he  based  on 
the  well  known  differential  equation  of  the  characteristics 
for  an  axially  symmetrical  isentropic  flow,  namely: 

t  d0  -  ~  cotya  +•  Siry*.  sin  0  -y-  -  0  (1) 

where  the  top  signs  refer  to  the  differentiation  along  the 
left  running  (L)  -  and  the  bottom  -  along  the  right  run¬ 
ning  characteristic  lines  (R  -  see  figure  2). 


Replacing  the  velocity  W  by  its  components  U  and  v, 
and  the  element  dl  by  the  expression  ±  > 
after  some  transformations,  one  gets: 
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dv+  cot  (yu+0)  + 


sia1^ _ 

Stg^+Q)  slg(yu-0) 


(1') 


In  order  to  simplify  the  above  equation,  two  assump¬ 
tions  will  be  made;  these  are  approximately  equivalent  to 
the  linear  theory  of  V.  Karman  and  Moore.  These  assumptions 
are: 

(i)  -  »  0  ,  which  means  that  angles  0  must 

be  relatively  small,  and  Mach  numbers  not  too 
great. 

(ii)  -  jj.  ~  ,  that  is,  the  Mach  numbers  are  al¬ 

most  constant  and  equal  approximately  to 
for  the  undisturbed  flow. 

The  assumption  (i)  implies,  obviously,  that: 

/X±  0  =!  yu. 

which  enables  to  present  equation  (1')  in  a  very  simple  fi 
nal  form: 


d(vy)  *  ft  y  dll  =  0  (2) 

The  assumption  (ii),  which  to  a  certain  point  is  the 
consequence  of  (i),  means  that  all  the  Mach  lines  may  be 
treated  as  mutually  parallel  straight  lines.  This  fact 
makes  the  integration  of  the  equation  (2)  considerably  eas, 

ier. 
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3.  THE  INTEGRATION  OF  THE  SAUER-HEINZ  EQUATION 

In  order  to  perform  the  integration  of  equation  (2), 
let  us  choose  a  certain  number  of  arbitrary  points  on  the 
body  surface  and  draw  therefrom  left  (L)  and  right  running 
(R)  Mach  lines,  parallel  to  the  Mach  lines  of  the  undis¬ 
turbed  flow.  This  procedure  will  create  a  "net"  of  Mach 
lines  which  we  shall  call  a  Mach  net. 


To  facilitate  the  calculation,  let  us  number  the  lines 
of  the  Mach  net  in  the  manner  shown  in  figure  3*  Points, 
which  are  the  intersections  of  the  left  Mach  line  with  num¬ 
ber  L  and  of  the  right'  one  with  number  R,  will  be  .designated 
oy  (L;R). 

Assuming  that  the  distances  between  neighboring  points 
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in  the  Mach  net  are  small  enough,  we  may  write: 

H«)  \ 

jy  dU  -  -jr  (yL-,n+  ")  ( U.L.*  -  - 

-  2^  (yi.;l».  +■  Vl;R-i)  (  U-  l.;ft  —  U.  l'R-1  ) 

and 

(tiR) 

J  ydU  cr  (  y R  +  - 


thua,  the  result  of  integrating  (2)  may  be  written 

(vy)c;*-(vy)u<M=  yMYsR  +  7 

(Vy)f,R~(V^)w-1-,R="^  t)/u-(;R')(u.LlR  -  Uvi;r)  (3  t  ) 

the  boundary  condition  is: 

vK;lt  =  (U.+  a)  tan  0K  2  (Uc*  +  u)  0K  v  1 


where  0K  denotes  the  slope  at  point  (K;K). 

The  problem  of  determining  flow  velocities  in  the  net 
points  will  be  now  discussed  in  four  parts  according  to 
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the  position  of  the  points, 

.4*  THE  DIFFUSER  NOSE  POINT 

Let  us  choose  two  points  on  the  same  Mach  line:  P  on 
the  diffuser  surface  and  Q,  in  the  undistruhed  region  (see 
figure  4) 


In  order  to  calculate  the  velooity  up  it  is  possible 
to  use  aquations  (5')  and  (J1*) 

Solving  these  equations  and  assuming  that  P-*-(0;0) 
and  ^-—(OjO),  one  gets  immediately: 

«■>  =  -  (Tferu«  “  -fu“  “*> 

where  the  subscript  n0"  refers  to  the  point  (0;0)»  The 
equation  (4)  is,  of  course,  equal  to  the  result  for  a  two 
dimensional,  flow,  deflected  through  an  angle  0O  •  Obvious- 
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ly,  this  must  have  been  expected,  considering  that  the  prob 
lem  in  the  close  neighborhood  of  the  nose  point  is  indeed 
two  dimensional. 

5.  POINTS  ON  THE  FIRST  MACH  LINE  (L=0) 

In  order  to  solve  the  problem  for  the  first  Mach  line, 
i.e.  the  Mach  line  with  L=0,  let  us  choose  on  it  two  points: 
(0;R)  and  (0$R-1)  (see  figure  5) 


On  the  right  Mach  line,  starting  from  (0,R),  let  us 
assume  a  point  CL  Recalling  that  the  flow  is  undisturbed 
at  Q.  and  assuming  that  Q.  is  approaching  (0,R)  (therefore 
y0.R  ),  one  gets  immediately  from  the  right  charac¬ 
teristic  equation: 

-  “  P  y  01P.  U.0;  R. 


(5) 
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It  is  evident  that  this  relation  is  valid  for  all 
points  on  the  first  Mach  line, and  particularly  for  (0,R-1), 
which  yields: 

—  — P>yo;R-l  ^OiR-l  (5») 


Substituting  values  (5)  and  (5‘)  into  the  equation  (3), 
one  can  calculate  the  value  of  u  0;R  as  a  function  of  u0;Rrl 


U.  0;R 


22&.±Syj!fcL.  U  ojR-i 

3  y»;R  4-  yo-,n.-i 


(6) 


It'  is  possible,  however,  to  obtain  the  value  of  u0.R 
as  a  function  of  y0.R  only.  For  this  purpose,  let  us  call 

Au  =  U.0;R  —  U0;  R-l 

^  y  ~  y°;p.  ~  y.,H 

The  equation  (6)  becomes  now 


u. 


o  -  R 


-  (u,,-Au) 


(7) 


Assuming  that  Ay  -*-0  and  neglecting  the  second  or¬ 
der  terms,  the  following  differential  equation  is  obtained 
(omitting  the  subscripts) 


d  u  _  1  dy 
u  “  2  y 


(7 ' ) 
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whose  solution  is: 

U.=  [X-1*U»  (8) 

.  y°. 

where  u0  and  y0  are  the  values  at  the  diffuser  nose. 

6.  POINTS  ON  TEE  DIPPUSER  SURFACE 

As  to  points  (K;K)  which  lie  on  the  diffuser  surface, 
admitting  that  the  flow  Is  reaching  them  with  the  right 
Mach  line  (see  figure  6)  equations  (31  )  and  (3'1)  must  be 


\  / 
\  / 


used.  After  solving  them  and  eliminating  (vy)*-!,*  by  means 
of  equation  (3),  one  gets 

—  Aw  Bk',K  +  Ck.K  Ole  (9) 

where  A,  B  and  C  are  certain  geometrical  expressions  depend 
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ing  on  positions  of  points  (K;K),  (K-1;K),  (K-l;K-l)  and  on 
slopes  of  the  diffuser  surface  at  points  (K:K)  and  (K-l;K-l)t 
Their  values  are: 


A  -  +  ^yK-<jK  +  y  > 

(l+y5)  Vk:k  +  y*-i;K 

(l-  iHy  w  +  yx-iix  I 

B  K/K  *"  .  ?cv\ 

(l+  pJVK-.K’  +•  Vk-IJK 

(n-fH* +y«,. 

Naturally,  assuming  a  Mach  net  such  that : 


(10) 


yo-.R-yojR-i  -  y0,M- yo,*-i  -  -  =  c°nst  =  Ay 


and  denoting:  yK.K  -  =  (Ay)K  the.  values  of 
A,  B  and  £  may  be  presented  in  a  form  more  convenient  for 
numerical  calculations,  i.e.: 


A«^ 

8k;k 


C 


K’iK 


Jn  k;k  +  2Ay  -3(Ay)K  j 

Z(l+  fl)yK-,K+Ay-(Ay)K 

1  (l~  ^ )  )/k)k + Ay  ~lji~  )(Ay)n  > 

Z(l4  |»)  yK.K  +  Ay  (Ay) k 

2(1+  !*)yK;K  +Ay  -  (Ay)  K 


(10'  ) 


As  to  the  point  (1;1),  using  the  more  accurate  equa¬ 


tion  (51  )  instaed  of  (3),.  one  gets: 


u 


v.i  - 


(U  4=) 


U 


O|0 


Ay  ~  (Ay)<< 


2(l+  p1)  y<«  +  Ay-  (Ay)< 


Uoji 


(9' ) 


7.  THE  INTERIOR  POINTS  OP  THE  NET 


Solving  equations  (3)  and  ( 3 1  )  for 


Uu*  = 


-  +  ^ L)*-<  4  2.  P  (>;»  V,*)  ;R 


|-P(lyi;R+  yt)R-1+  y l-1;R  ) 


(11) 


writing 

(^yJu-r.R-  (vy)MM  =  [(vy)L-1;ft-  (vv)u-1;r.i]  - 

and  utilizing  directly  equations  (3)  and  (3')  to  calculate 
the  above  expressions  as  a  function  of  u,  upon  eliminating 
the  values  of  (vy)  we  get  finally: 
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U-lr  =  D  L;R  U-L-V,R.  ~  b  t--,R 


where  D,  E,  and  F  are  certain  geometrical  expressions,  de¬ 
pending  on  the  Mach  net  assumed.  Their  values  are: 


3  V  kiiS  ±_¥kliS=l 
2  >l;r  +  +  y^i;R 


3  v-L-.R-i  ~t~  Vl-i’ik 

lyujR.— yu-,RH  +yL-<;R 


> 


(12) 


r  -  3  yw,R-i+3yt-ii»  -Zy.ua. 

2.yu-,i»+  yt.-,R-i  •hyL'-ijR. 


Assuming  the  same  Mach  net  as  in  the  former  section, 
the  values  of  D,  E,  F  may  be  presented  in  a  simpler  form: 


n  _  ^1-yi.iR  +ZAy  ~  3  (Ay)  u 
ULiH  '  4yu,,.  -  (Ay)  l 


-  4vm,*  -  2  Ay  -  (Ay) l 

4  y  hr  -  ( A  y)  u 


-  4y mr  ~  3  (Ay )  l 
4yf,R  -(Ay)L 


(12 ») 


Calculating  particularly  the  points  of  the  ''second" 
Mach  line  (i.e.  with  L  =  1),  in.  order  to  eliminate  the 
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values  (vy),  it  is  possible  to  use  equations  (5)  and  (5') 
-more  exact-  instead  of  (3).  Thus: 

(Vi)o;R  -  (Vy)l;R-1  =  j3yo;R-1  U-0;R-1  -  Ji  yo;R  ^O-.R  + 

+■  Y  P(yi;R-i  +  yo',R-i)(^  bR-1  ~u-o;R-<i') 

Substituting  this  expression  into  equation  (11),  and 
solving  it  for  ui;R  ,  one  gets: 

b*-1;R  =  E-t;R  WV,R-1  "*■  (Uo.g^-  U-O-r)  (ll") 

where  the  value  of  E  1;K  is  determined  by  formula  (12')  and 


G 


»\R 


_ Vo'.»  ~  y^.R 

2y<;R+  yi-,R-i  +yo;R 


Ay  -(Ay)< 
4y«;*-(Ay)< 


(12"  ) 


As  to  expression  -  u0;r)  in  equation  (11"), 

it  is  possible  to  present  it  in  the  form  -see  equation 
(7)-: 


u°iiM  ~  n-oiR  —  —  U. 

^■yois-Ay 

The  equations  (4)  and  (8)  enable  to  calculate  the 
u-velocity  directly  from  the  geometrical  positions  of 
points,  and  the  equations  (9),  (9' ),  (11*),  (11")  -  from 
the  known  data  at  former  points  of  the  Mach  net. 
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8.  A  PRACTICAL  PROCEDURE  TO  CALCULATE  THE  VELOCITY  DISTRIBU 
TION  ALONG  THE  SURFACE  OF  A  DIFFUSER 

The  practical  procedure  may  be  presented  as  follows: 

On  the  first  Mach  line  of  a  diffuser  one  chooses  cer¬ 
tain  number  of  equidistant  points  (including  the  nose  point) 
so  that  the  Mach  net,  starting  from  these  (see  figure  3) 
covers  the  entire  region  where  the  velocity  distribution  is 
to  be  calculated. 

After  finding  from  equations  (10*  )  the  coeficients  A, 
B_  and  C  for  all  surface  points,  and  D,  E,  F  from  equations 
(121  )  for  all  interior  points,  one  finds  the  velocities  u 
at  different  net  points,  in  the  following  order: 

(i)  -  nose  point  -  equation  (i|) 

(ii)  -  points  on  the  first  Mach  line  -  equation  (8) 

(iii)  -  point  (l;l)  next  to  the  nose  on  diffuser  sur¬ 

face  -  equation  (9' ) 

(iv)  -  points  on  the  Mach  line  L  =  1,  issuing  down¬ 

stream,  from  point  (1;1)  -  equation  ( 11” ) 

The  procedure  described  in  items  (iii)  and  (iv)  is 
next  repeated  in  the  same  order  with  respect  to  other  sur- 
face  points  -  equation  (9)  -  and  lines  L  =  const  -  equation 
(11')  -  issuing  therefrom,  till  the  last  net  point  in  the 
diffuser  surface  is  reached. 

It  is  necessary  to  mention  here  that  in  order  to  ob¬ 
tain  the  velocities  u  for  n  points  on  the  diffuser  sur- 
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face,  they  must  be  determined  also  for  ^h(n-l)  remaining 
net  points.  So,  on  the  one  hand,  the  more  dense  is  the 
Mach  net,  the  more  accurate  the  results  will  be,  but,  on 
the  other  hand,  the  amount  of  work  required  increases  con¬ 
siderably. 
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Part  two  -  STRAIGHT  LINE  DIFFUSERS 

9.  ADAPTATION  OF  THE  METHOD  TO  STRAIGHT  LINE  DIFFUSERS. 

In  case  of  a  straight  line  diffuser  the  amount  of 
work  may  be  reduced,  assuming  that  the  distance  between 
two  neighboring  net  points  on  the  surface  of  the  diffuser 
is  constant  (see  figure  J).  Thus,  the  x  co-ordinates  of 
these  points  will  be  K.m,  where  K  =  0$1;2;  and  m  is 
certain  arbitrary  length  on  which  depends  the  net  density 
and,  therefore,  the  accuracy  of  the  solution. 


It  results  from  geometrical  relations,  that  (-»-  sec¬ 
tion  U) : 
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y™  =y.-^d-P«)i-  + 

(Ay)K  -  const  —  mtan  0  —  m.0 

Ay  =^(l+po) 

Substituting  these  values  into  equations  (8),  (9), 

(9'  ),  (11'),  and  (11*),  and  dividing  them  by  u.0  one  gets: 
0  For  the  first  Mach  line: 


* 


®  For  the  dl  if  fuser  surface: 


4.  ±-m 
a«  "  (1+f )Vks«+tjf 


m. 


tUim  _  frfWL  1^+^ITl  UlMiiH  i 

u‘  (n-f)>:*+ Uo 


.  m9  (14) 


particularly,  for  K®1,  one  may  use  the  more  exact  equation: 


Up-,  A 

Uo 


d4') 
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(c]  For  the  interior  points: 


and  particularly  for  L  =  1,  one  may  use  more  exact  equa¬ 
tion: 

Ui;n 

U0 


where  the  coefficients  of  velocities  ( u/u.0)  ,  in  spite  of 
their  apparent  complication >can  be  determined  numerically 
for  a  given  diffuser  and  a  given  flow  quite  easily  and  rajD 
idly. 

Three  graphs  are  enclosed  -  1  to  3  -  (see  -pg.60/6?.)  , 
which  show  the  u  -velocity  distribution,  found  by  means  of 
the  above  described  method.  The  angles  0  are:  0,025  radi¬ 
ans,  0,05  radians  and  0,125  radians;  is  assumed  equal  to 
1,  and  the  segment  m  varies  from  0,2  y„  for  initial  points 
to  0,8  \/v  for  last  points  of  net.  The  computations  were 
performed  in  the  region  0  ^  x  ^  "16  y0  .  The  dependent  va¬ 
riable  is  u/u? 


22 


10.  AN  OUTLINE  OF  THE  ANALYTICAL  SOLUTION  IN  THE  CASE  OF  A 
STRAIGHT  LINE  DIFFUSER. 

The  method  described  in  the  preceding  section  fur¬ 
nishes  a  possibility  of  finding  in  a  simple  way  numerical 
results  for  the  velocity  distribution  along  the  surface  of 
a  straight  line  diffuser.  It  doesn't,  however,  determine 
the  analytical  form  of  this  solution,  nor  its  relation 
with  the  angle  9 . 

In  order  to  answer  this  question,  let  us  assume,  for 
the  sake  of  argument,  that  the  diffuser  extends  downstream 
to  infinity.  Obviously,  in  the  case  of  supersonic  flow, 
this  does  not  restrict  the  generality  of  a  solution. 

Let  us  call  a (x)  a  function  of  x  which  represents 
the  u -velocity  distribution  under  investigation  along  the 
surface  of  the  diffuser.  The  value  of  this  function  at  the 
nose  point  and  at  infinity  can  be  immediately  established. 

In  fact,  at  the  nose  point,  as  it  was  already  shown 
in  the  preceding  section,  the  velocity  u(o)  (which  will 
be,  as  earlier,  indicated  by  u0  )  is  equal  to  the  corre¬ 
sponding  velocity  in  a  two  dimensional  deflection,  that  is: 

|  Ucc 

As  to  the  velocity  u.  (<»)  (which  will  be  indicated  as 
Uoo  ,  it  is  possible  to  use  a  more  intuitive  similarity  or 
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"scale  effect"  reasoning.  Consider  at  first  that  the 
length  of  the  diffuser  is  limited  and  equal  to  a  (see  fig¬ 
ure  8).  The  radius  of  its  nose  point  is  y0  .  If  we  now  ex¬ 
tend  a  to  infinity,  but  choose  the  scale  of  the  drawing 
such  that  this  length  remains  always  the  same,  then  y0  (in 
the  figure)  will  tend  to  zero,  becoming  zero  when  a  =  oo . 

In  this  case  we  obtain  a  cone.  Thus  the  velocity  on  the 
surface  at  infinity  for  a  given,  supersonic  flow  is  unique¬ 
ly  defined,  since  the  velocity  of  the  surface  of  a  cone 
does  not  depend  on  the  position  of  the  point  but  only  on 
the  angle  Q . 


One  can  obtain  this  velocity  from  appropriate  plots, 


or  from  cone  tables,  or  approximately  (for  small  values  of 
0  ): 


In  order  to  establish  the  form  of  the  function  u.(x) 


for  values  0<x<oo  ,  it  is  convenient  to  assume  an  euxi- 
liary  function: 


g(x)  -  Hop 
Uo  —  g  oo 


(16) 


which,  independently  of  9  and  J3  ,  for  x  =  o  and  x  =  <» 
admits  the  following  values: 


°*(o,e,p) 


U.0  ~  U.  QQ 

go  -  u-oo 


1 


—  g«>~  u.ot> 

U-o  —  Uoo 


0 


The  behavior  of  the  function  (X  for  other  values  of 
x  is  7  for  the  time  being  -  unknown.  It  will  be  invest igat 
ed,  however,  in  the  next  sections.  Once  the  function 
o'  (*.e>£)  is  established,  finding  of  the  function  u(x') 
becomes  trivial,  because  from  equation  (16): 

g(x)  -  Uoo  +■  (Uo  -  U00)-O'(x,9,(3)  (17) 

We  see,  from  equation  (17),  that  the  function  o'  has 
a  very  simple  physical  interpretation:  it  indicates  the  man 
ner  in  which  the  transition  from  two  dimensional  flow  at 
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the  nose  point  (  <r=  1  u^)=  U.0  )  to  three  dimensional 

conical  flow  at  infinity  (  O' =  0  u.(xA  =  U-oo-)  is  per¬ 
formed  along  the  surface  of  straight  line  diffusers, 

11.  DETERMINATION  OF  THE  LIMIT  OF  THE  FUNCTION  O'  AS  THE 
ANGLE  0  APPROACHES  ZERO. 

When  0  =  0,  the  flow  is  undistrubed  and  the  func¬ 
tion  or  has,  in  this  case,  an  undetermined  value: 

Cfc.o.jO  =  § 

Regardless  of  this,  as  will  be  proved,  the  limit  of 
this  expression  exists.  Let  us  indicate  it: 


To  prove  its  existence,  let  us  note,  that 


U_m  crCx.Q.p) 


Um  -  Uo*~  = 

Q-—0  —  U-oo 


Urn  7=^  -  lim 

9-*0  Up _ Q-*»  Q 


LL  QQ 

H7 


1 


lim 

9-0 


Uqd 

U.0 


Using  L’Hospital  rule,  we  have! 


26 


(im 

e-o 


L  im 
a— o 


4 


pe 


lira  (pe)  =  0 


from  which 


U_m  c rM^) 


l  im. 

e-*o 


<±iA 

\Xq 


(18) 


The  function  ul>0/u.o  can  be  calculated  numerical¬ 
ly  by  means  of  equations  (13)  5  (lit),  ilk'),  (15)  and  (15'). 
We  recall  that  all  these  equations  are  determined  when  0-~O 


and  have  the  following  limits  when  0=0  : 

[A]  For  L  =  0: 

tel  =  [i+  W* 

(19) 

1] 

For  -L  =  R  (  =  K 

tO)  : 

fikAUS.' 
L  u,  . 

-  7  [jAiriUsl  f  U.K-tiK-<l 
’  L  L  U0  J  L  U0  J 

(20) 

and  particularly  for  K  =  1 

fe*-] 

a  m,* 

I  if  [U7J 

!+?•* 

(20') 

0  For  L*Rj  L*0 
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U-l-,r 

I  Z  +  (R-L-+4)m* 

Ul-i-.r 

1 

? 

T 

3 

*. 

UL;R< 

M 

U.L-1-.R-1 

U.Q  _ 

Z-t-(R-L)rn* 

_  „ 

2  +  (R~L)m*l 

L  U0  J 

LU  . 

(21) 


and  particularly  for  L  =  1: 


ru4;Ri 

|_  2+(R-l)m*  i 

[U-t-.R-f 

1 1  £ 171  *  f 

fUoiwl 

[Wl\ 

L-a0  J 

2  +  (R-i)m* 

1  J 

Z+(R-l')  m*  1 

L  u0  J 

L  Ij 

(21*) 


where: 


Since  all  these  equations  are  determined  it  is  evi¬ 
dent  that  the  limit  (18)  exists  and  may  be  numerically  cal¬ 
culated  at  point's  Km  by  solving  for  u-k;k  /  u0 
Moreover,  it  is  clear,  from  the  structure  of  the  above  equa 
tions  that  this  solution  is  a  function  of  a  parameter  m*  = 
=  m/[iyD  .  Thus,  introducing  a  new  variable-: 


to  each  point  x  =  Km  on  the  x  axis,  there  will  correspond  a 
point: 


5 


K  m* 


on  the  ^  axis.  The  values  u*.*/  u0  will  be  now  subordi¬ 
nated  to  the  points  Km*  but  as  they  depend  only  on  the 
choice  of  m*  the  function  o'  can  be  expressed  as  a  func¬ 
tion  of  only  one  variable  \  Thus: 
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Um  CT  (x.Q.p)  =  |(0  =  {(%,,) 


12.  NUMERICAL  EVALUATION  OP  THE  FUNCTION  O' [\) 

Substituting  into  equations  (19-21)  any  numerical  val 
ue  for  m*  and  solving  for  uk,k//u.,>  we  shall  obtain  the 
values  of  the  function  CT(0  at  the  points  \  -  K  m*  ,  The 
following  tables  give  the  values  of  cr  as  a  function  of  \ 
for  m*  =  0,V  and  for  m*  =  0.2  in  the  range  0  k  \  L  4 


1 

<n*) 

0 

i 

1 

0,2 

i  -  , 

- 

0,h 

0.825 

0,6 

0.756 

- 

0,8 

0,69k 

0.693 

1,0 

0..6U0 

- 

1,2 

0.592 

0.592 

1,1; 

0,5k9 

- 

|Q 

0.511 

0.511 

1,8 

0.U77 

2 

0,kk7 

o.kbl 

1 

m*  =  0,V 

o.l;l9 

- 

2,1; 

0.391; 

0.391; 

2,6 

0.371 

- 

2,8 

0.351 

o.35l 

E 

0.332 

- 

3,2 

0.315 

0.315 

3,1; 

0.299 

- 

3,6 

0.285 

0.285 

3,8 

0.271 

- 

k 

0.259 

0.259 
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Substituting  for  m#  values  less  than  0,2  one  will  not 
obtain  any  changes  in  the  first  three  decimal  points  of  the 
function  CT(^)  *  The  result  obtained  for  m*  =  0,2.  can  be 
therefore  considered  as  a  final  one  for  the  first  three  de_c 
imal  places. 

For  greater  values  of  $,(^>4-),  it  is  necessary  to  in¬ 
crease  the  value  of  m*  .  Although  such  results  for  u  K.  K 
at  small  values  of  _K  are  always  determined  with  relatively 
greater  error,  as  K  increases  a  tendency  towards  '‘self  cor¬ 
rection"  is  observed.  Thus,  combining  different  results 
(with  different  values  of  m*),  it  is  possible  to  reproduce 
relatively  exact  values  of  cr(^). 

These  values  are: 


1 

cnj) 

5 

0.208 

6 

0.175 

8 

0.131 

\ 

<y<.\) 

10 

0.104 

20 

0.051 

40 

0.025 

From  the  last  results  a  tendency  can  be  distinctly 
seen  for  the  function  CT(^)  to  approach  \~i  for  great  val¬ 
ues  of  l  , 

Graph  4  showfe  the  shape  of  the  function  CT  (\)  in  the 
region  0  ^  \  A  16  . 
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13.  THE  SUCCESSIVE  ANALYTICAL  APPROXIMATIONS  OF  THE  FUNC¬ 
TION 

Let  us  solve  the  system  of  equations  (19-H11)  for 
UK;K  (K  =  l;2.;3;.„)  ,  without  substituting  any  numerical  val. 
ue  for  m*  .  These  solutions  will  be,  thus,  functions  of 
m*  ;  let  us  indicate  them  T*  (m»)  ;  Tj ^*)  ;  ...  .Generally, 
each  function  VK (r^*)  determines  approximately  (with  accu¬ 
racy  increasing  with  K  )  the  val-qe  of  the  function  C(^)  at 
point  l  =  K  m*  .  Thus,  for  each  concrete  value  of  K  we 
have: 


But,  as  the  value  Km*  (for  concrete  K  )  is  not  fixed  and 
may  admit  any  arbitrary  value  )  ,  we  have: 

-  MI/k) 


So,  the  sequence  of  the  functions:  ^  ;  (_ | /2.)  ; 

Tj  ...  represents  the  successive,  more  and  more  accu 

rate  approximations  of  the  function  <r($) 

In  the  case  of  K-i  one  obtains: 


Til*)  - 


zil  +Vi+-H 
(i+tOVT+F 
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Expanding  the  value  of  \f[  +  \ $  1  into  power  series 
and  keeping  only  the  first  two  terms,  one  obtains  immediate 

lys 

°h>  ~  Tl(°  ~  (TTiT7  (22) 

as  a  first  approximate  of  the  function  O' (O' 

Choosing  K=  1  ,  the  functions  T  become  at  once 
very  complicated.  For  instance,  already  for  K.  =  2.  ,  there 

is: 


\3-  4^-+  32^+256 

256(l+^)M1+i^r 


(25) 


For  practical  purposes,  the  approximation  given  by 
equation  (22)  is  quite  adequate.  Being  very  simple,  it  pro 
vides  at  the  same  time  quite  -good  accuracy,  especially  in 
the  range  0  4  \  4  h  ;  this  may  be  demonstrated  by  a  com¬ 
parison  with  the  yalues  of  the  function  cf(^)  given  in 


section  12 
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$ 

cr<p 

(1+K)"Z 

0 

i 

1 

0,2 

0.907 

0.907 

o,4 

0.826 

0.826 

0,6 

0.756 

0.756 

0,8 

0.69k 

0.694 

l 

0.614,0 

0.640 

1,2 

0.592 

0.592 

1,4 

0.549 

0.549 

1,6 

0.511 

0.510 

1,8 

0.477 

0.476 

2 

0.447 

0.444 

o'(S) 

2,2 

0.419 

0.416 

2,4 

0.394 

o.39l 

2,6 

0.371 

0.367 

2,8 

0.351 

0.346 

3 

0.332 

0.327 

3,2 

0.315 

0.309 

3,4 

0.299 

0,292 

3,6 

0.285 

0.277 

3,8 

0.271 

0.263 

4 

0.259 

0.250 

The  behavior  of  the  function  cr(^).for  large  values 
of  l  may  be  investigated  by  means  of  its  higher  approxi¬ 
mations:  31(1/2)  5  3"-j  (1/3 )  ;  ....  J  it  can  be  shown  that  they 
all  approach,  (for  large  values  of  ^  ),  the  functions  a ^ 

where  a.  varies,  depending  on  K  ;  as  to  the  function  3  (\) 
itself,  it  seems  from  numerical  values  given  in  section  12, 
that  a=  1  ,  i.e.  that  <T(^)  approaches,  for  large  values  of 
\  the  function  V1  (  this  will  be  confirmed  in  next  sec¬ 
tions). 

The  comparison  between  the  function  cr(^)  and  its 
first  approximation  [l+|\]  1  for  large  values  of  \  ,  shows 
that  both  functions  approach  zero;  however,  the  former  does 
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it  like  the  function  $-1  while  the  latter  as  16  |  1  ..Thus, 
the  absolute  error: 

(l+iO'* 

which,  as  can  be  established,  is  always  less  then  0.025, 
will  also  approach  zero  for  large  values  of  ^  .  The  rela¬ 
tive  error  will,  of  course,  tend  to  infinity.  For  practi¬ 
cal  computations  of  the  u. -velocity,  however,  with  certain 
restriction,  it  will  have  no  major  effect  (see  section  17). 
Graph  5  shows  the  comparison  between  the  functions 
+  and  ^,_1  in  the  range  0  4  {  4  16 

1 U*  THE  OPERATIONAL  FORM  OF  THE  FUNCTION  CT(^) 

It  is  possible  to  obtain  the  exact  analytical  repre¬ 
sentation  of  the  function  O' [\)  by  means  of  operational 
calculus. 

For  this  purpose  one  will  use  the  transformations 

CO 

F (p)  =  P  J fix)  ^dx  (210 

0 

which  will  be  denoted  by  the  symbols: 

F(p)  -  P[j(x)] 

The  inverse  transf ormatibn  will  be  denoted  as: 


3h 


f(*)  -  P  [F(p)] 

As  it,  is  well  taown,  the  inverse  transformation  may¬ 
be  determined, under  certain  restrictive  conditions, by  means 
of  the  formula.: 


P"[M 


Xo  +  lOO 


d* 


Xq-IOO 


(25) 


where  x0  must  be  chosen  so  that  the  line  Re  (i)  -  Xo  Pas“ 
ses  to  the  right  of  all  singularities  of  the  integrand. 

The  function  o*(|)  may  be  determined  by  solving  the 
differential  equation  for  the  perturbation  potential  in 
three  dimensional  axially  symmetrical  flow: 

+  rO  (26) 

by7  ^  y  i> y  I3  J)xl  w 

with  boundary  conditions 

(V\  =  00  -  $y\sco  ~  0 

C^V^ex  -  (|^)v=y.+ex  =  9U°°  J 

(once  it  is  assumed  Q-*-0  ,  the  second  term  of  the  sum 
y0+9x  will  be  omitted). 
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Transforming  equation  (26)  by  means  of  equation  (22+) 
with  respect  to  the  variable  x  ,  and  taking  into  consider, 
at ion  that: 

-P  P|/]  ~  ^PM 
■P  B  If]  =  757pM 

p  Hi]  =  P!{  PlP  M  -  p‘t„+o  -  ?($•),<=*„} 

and  that  the  choice  of  the  co-ordinate  system  is  such  that 
for  x  =  -i-O: 


<f  =  0 


%*-  =  u.  =  0 
ox 


one  obtains  the  equation: 


^P[<f]+7^PW-pVP[<f]  =0  (28) 


the  solution  of  which  is: 


P  [(p]  -  Cilo(fJyp')  +  C*  K0  (J3y  f>) 


(290 


Submitting  the  boundary  conditions  to  the  transforma¬ 
tion  ?,  one  obtains: 

%]),...=  0 

=  P[0U*]  =  9U“ 

Applying  the  above  conditions  and  taking  into  cons id 
eration  that 

^PM  =  C1ppI1(pyp)-ClppK1(pyp) 


one  gets: 


Ci  =  0 


c2  = 


,  8  Uoo 

PP  Kt  (jiy.p) 


and  finally: 


PH 


_ K.°  (pyp) 

PpKi(py„  p) 


euoo 


(30) 


The  operational  representation  of  the  velocity  u  (x,y) 
with  respect  to  x  ,  is: 
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y  [u(^v)]  -  P  [|f ]  -  p  P  Pf]  - 


Ko(B  V  p)  Q  Hop  ._  KoCftyp')  ^ 
Kl(^-y0p)  P  Kl(py0p) 


(31) 


The  last  result  enables  us  immediately  to  establish 
the  o;  erabional  form  of  the  function  O' ({)  .  In  fact,  com¬ 
bining  equations  (18)  and  (31) ,  one  gets: 


P  bix.pd  =  p  [•“£?*>]  = 

i.e«,  the  final  result. 


(32) 


15 o  CALCULATION  OF  THE  FUNCTION  CT(Xjp)  FROM.  ITS  OPERA¬ 
TIONAL  FORM. 

In  order  to  find  the  function  o'U.jV)  one  has  to 
perform  the  operation: 

P'1 

With  this  purpose  the  operation  will  be  performed 


Kolft-ypp) 


Kt(Ps0p)J 


(33) 


first 
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i.e.,  the  integral  will  be  calculated 


(33’) 


X-+  Coo 


2JTi 


KoM  e 


xz 


Kite)  i 


dz 


X.-  too 


(33") 


and  next,  "shift"  formula  of  operational  calculus  will  be 
applied. 

As  it  is  known,  the  functions  K0(2.)  and  Ki(i)  have 
a  logarithmic  branch  point  at  x  =  0  ;  besides  this  they 

have  no  other  singularities.  Moreover, ''the  function  Ki  (z) 
does  not  have  any  zeros  in  the  region  |argH|<rc  so  the  in¬ 
tegrand  of  (33")  is  holomorphic  therein. 

As  a  result,  we  may  write: 


where  1C  =  ABCDEFG-HA  is  the  path  of  integration,  shown  in 
figure  9* 
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'  <0  —  0 
r  —  0 

< 

£  — -  0 
k  R  — oo 


then  the  path'X  may  be  divided  (in  the  limit)  in  the  fol™ 
lowing  manner? 

3C  =  JC+  O  C'  +  A  +  ff 


where: 


ko 


£  -  is  the  path  on  the  y  axis  from  -  co  to  +  oo  leajr 

ing  the  point  z  =  0  to  the  left  by  means  of  the 
semicircle  with  r— 0 

C'  -  is  the  semicircle  of  radius  R  —  oo  traced  from  the 
point  z  =  0  to  the  left  of  the  y  -  axis 
C"  -  is  a  circle  of  the  radius  r— 0  traced  from  the 
point  z  =  0 

M.  -  is  the  path  from  -  oo  to  0  on  the  x  axis,  con¬ 
sidering  that:  -x  =  reLir 

-  is  the  path  from  0  to  -  oo  on  the  x  axis,  con¬ 
sidering  that:  -x  =  re"lTr 
The  behavior  of  the  integral  (33)  along  all  these 
paths  will  now  be  investigated, 
a)  We  have: 


r  Ko  (z) 

J  K1(Z) 

c' 


(35) 


which  may  be  proved  by  means  of  Jordan’s  theorem,  taking 
into  consideration  that: 


Urn 

l -—CO 


Ko(z) 

zKife) 


l  tm 


-  0 


b)  The  integral 


U 1 


because 


J 

C" 


Kife) 


XH 


dz  =  0 


f 

c* 


Kz e*z  dz 
Kito  * 

Urn 

?-»>o 

-  2jrerx  Um 

H-*0 

*=- 

1 

c 

1 

t  0 

12  |  1 

7  +  T  laT 

Koto 


Kito 


c)  In  order  to  calculate  the  Integrals 


KoM 

Ki  to  ?. 


Koto 

Ktto 


M 


dz 


one  uses  the  following  relations-: 


K  o  feetm)  =  K0 (to  ^  151 1 o to 
Mze^’1)  =  -<iU)  +  i-Kj.to 


(36) 


therefore,  after  a  simple  transformation: 


i)2 


J 


KoMUd)  t  KjjtUoM 
K/(t)  +  JT2  Ii(t) 


Taking  into  consideration  that 


Ko(t)  Ii(t)  +  Io(t)  “  ^ 


it  gives: 


w 

]  =  -2iuJ 


,  -tx 


dt 


Ktl(t)  •+*  Ttz  I^(t)  t* 


(37) 


d)  Finally,  as 


fj<aM  ex*d 

u  1  P'1 

r  Ko(P>i 

J  Kite)  i 

ZlTl  ~ 

L  J 

jC 


(33) 


considering  equations  (3 U  ♦  38),  one  gets: 

niELl  -  F  ,  e'*" _  dt 

Ki(p)J  J  lC(t)  4-  n1  V(t)  f- 

0 

After  applying  the  scale  effect  of  operational  cal~ 


cuius,  the  final  result  is: 


U3 


-  Ki(Pvp)- 


I 


tx 

e 

Ka2(t)  +JtlI4,Ct) 


dt 

t*- 


09) 


It  is  interesting  to  mention  here  that  for  y0  =  0 
the  value  of  e  U//^°  is  o,  and,  consequently,  one  ob¬ 
tains  0'(xR)=  0  ,  and  therefore  U(^  =  uoe+Or()<ip')(uo  +  U.oo)  =  U«>  5 

this  osult  is  obvious,  taking  into  considerations  that  for 
yQ  =  0  the  perturbation  body  is  a  cone,  and  that  is 

(  — —  section  10)  the.  u. -velocity  at  a  cone  surface. 

For  the  same,  reason,  when  the  curvature  is  infinite, 
one  has  two  dimensional  flow.  Therefore,  substituting  into 
equation  (39)  the  value  y 0  =  00  ,  one  must  obtain: 

OO 

( _ 1 _ <a  - 1 

0 

which  may  be  interesting  even  from  the  purely  mathematical 
viewpoint. 

16.  THE  APPROXIMATE  SOLUTION  FOR  O' (x,Jj)  FROM  ITS  OPERA¬ 
TIONAL  FORM. 

Considering  the  asymptotic  expansions  of  the  func¬ 
tions  Ko  and  K*  ,  one  obtains  for  large  values  of  t  : 


bk 


Koto 

Ki(z) 


63  + 

128 


Remembering  that 


P'W"] 


a  x 

n! 


the  expansion  of  the  function  C (x,p)  for  small  values  of 
x  (after  considering  the  effect  of  scale)  is: 


±\ 

46  L 


X 


i\. 


n_r^i4_ 

1024  L  p'yoi 


(4  0) 


It  is  now  possible  to  compare  the  present  result  with 
the  approximate  solution  given  in  section  13?  formula  (22). 
Expanding  this  expression  into  power  series,  it  results  in: 


(1  +  4[2Vo) 


A  _  1  _X_ 

1  2- (3* 


-W 

i6  Lpy0  J 


+ 


10 

1024 


-»  l 

—  1 


(40') 


Both  series  (40)  and  (40’)  diverge  only  from 
fifth  term. 


the 


for  small  values 


'~1  rKP(py0p)~ 

K-i(jiy0p) . 

of  p  one  obtains  the  asymptotic  approximation  of  the  func- 


Calculating  now  P 
;yn 

tion  CT^p)  for  large  x.  The  result  is: 


k5 


fiVo 

X 


oa) 


or; 


<T(t)  =  |  0*') 

which  was  already  supposed  in  section  13. 

17.  THE  APPROXIMATE  ANALYTICAL  SOLUTION  FOR  THE  FUNCTION 

u(x)  • 

As  known  from  section  10,  equation  (17),  the  func¬ 
tion  u  (x)  ,  which  represents  the  distribution  of  the 

U  -velocity  along  the  surface  of  a  straight  line  diffuser, 
may  be  presented  in  the  form: 

U0O=  •  (u0-  Uoo) 

where  the  symbol  u0  denotes  the  perturbation  velocity  at 
the  nose  point  of  a  diffuser  (equal  to  the  two  dimensional 
perturbation  velocity  deflected  through  the  angle  0  )  and 
Uqj  -  the  corresponding  conical  perturbation  velocity. 

For  small  values  of  angle  0  ,  remembering  that: 


hS 


it  is  possible  to  replace  o^x.e.p)  by  the  function  <T(X)p) 
In  this  case,  as  a  first  approximation  of  the  function  u.(x) 
one  obtains: 


U(x)  *  Uoo  +  0^(Xip)  (u0-  u„)  0j2) 

Though,  as  it  was  mentioned,  this  formula  is  valid 
first  of  all  for  the  values  9—0  ,  it  may  be  used  in  a 

relatively  wide  range  of  angles.  This  is  due  to  the  fact 
that  the  value  of  (Un-u^)  in  relation  to  u«,  decreases 
with  increasing  values  -of  9  . 

For  instance,  for  9  =  10°  (and  =  1): 


U  o  -  U  qq  _  a  / 

Uoo  ~ 

while  for  0  =  20  0  (  p  =  1): 

!*■.?  ~  u9b  —  o  6 
Uoo 

Thus,  when  the  angle  9  increases,  the  errors  due  to 
replacement  of  by  are  compensated,  be 

cause  &lso  in  view  of  the  fact  that  -  1  )  the  sec, 

ond  term  on  the  right  hand  side  of  equation  (Zjij.)  -  where 
the  errors  are  generated  -  is  relatively  much  smaller  than 
the  first  one. 

Finally,  replacing  O' ^  ^  by  its  approximation 


hi 


[1+-JM]"2  -  in  the  range  of  0  4  |  ^  Zj.  the  error  due  to  this 
is  absolutely  neglegible..  In  the  range  \  >  U  this  error  in¬ 
creases,  but  since: 


&(\)  ~ 


<  0,015 


(see  section  13),  it  may  influence  the  numerical  calcula¬ 
tions  only  if: 

(i)  the  diffuser  is  very  long,  and 
(ii)  the  expression  ~U:°u^<>g  is  very  great  (which, as 
was  said  in  this  section,  takes  place  only  if  the 
angle  0  is  extremely  small. 

In  conclusion,  remembering  that  \  =  ,  and  un¬ 

der  the  condition  that  (i)  and  (ii)  does  not  hold,  the 
first  approximation  for  the  distribution  of  the  u.  velocity 
along  a  straight  line  diffuser  is: 


^■(x)  —  *1  oo 


U.o  —  U.00 


(b3) 


In  the  case  in  which  0  is  very  small  (for  instance 
when  9  <  0.025  r.adians)  and  at  the  same  time  \  is  very 
large,  we  can  represent  by: 


ii8 

-  it  i  U.o  ~  Uot> 

-  Uoo  t  .  ,  m  the  range 

J£4,  i.e. 

0i3’) 

-  Uco  +  -^-TU-PP  Q  .  for  4;  i.e. 

The  enclosed  graphs  (6,7,8)  compare  the  results  for 
u(x)  obtained  by  means  of  equations  (13-15)  for  angles 
0,  =  0.025  radians,  0*.  =  0.05  radians  and  e5  =  0.125  ra¬ 
dians  with  those  obtained  by  means  of  equation  (Z;3)«. 

As  can  be  seen,  the  agreement  for  all  three  angles  is 
quite  good.  Using  the  method  of  characteristics,  the  compu¬ 
tation  of  takes  several  hours,  while  using  the  equa¬ 

tion  (U3) 1  this  time  is  restricted  only  to  minutes. 

18.  THE  DISTRIBUTION  OF  PERTURBATION  VELOCITY  IN  THE  SPACE 
AROUND  A  STRAIGHT  LINE  DIFFUSER. 

Taking  into  consideration  the  result  obtained  in  sec¬ 
tion  Ik,  equation  (31),  the  u.  velocity  around  the  straight 
line  diffuser  in  the  case  9-»0  may  be  presented  in  the 
following  operational  form; 


U9 


After  considerations  similar  to  those  of  section  15 
it  is  possible  to  obtain; 


00  _  tx 

u  (XiV )  _  f  Kolvt)  Ii(lA  +  Ki(t')  joCyt)  6-  ^J1  (U5) 

u «  J  Kd(t')  +  rt1 L.l(t)  t 


where  y  =  Y/yo 

Expanding  the  function  into  asymptotic  se- 

ries,  one  gets; 


Kq(ap) 

Ki(p) 


4  -(“-l)p  r 

a  e 


1 


1^3(x  -1  9  +  6  a  +  3  5  a*  D-i 

8a  ^  UB  a1  “ 


_  75  +  27a  +  55a*  +  249  a’ 

102.4  a' 


+ 


0|6) 


Submitting  each  term  of  (Z4.6 )  to  the  transformation 
P  1  ,  taking  into  consideration  the  "shift  effect"  and  sep¬ 
arately  the  "scale  effect",  one  finally  obtains; 


O'.O 


33  +  6v~1+9v~z  MM  2 

25  6  LjiyJ 


83+llv-1  +  9y-I4-2Sy  ~3  1 3  \ 

2.048  Ljiy*l  —J 


(U7) 
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where  x.=  x- j3 (y-y0) 

As  can  be  seen  from  figure  10,  the  variable  x  has  a 
very  simple  geometrical  meaning.  It  represents,  namely, the 
distance,  on  the  horizontal  axis,  between  the  point  a  (x,y) 
and  the  first  Mach  line.  It  is  possible  therefore,  to  con¬ 
sider  it  as  the  shifted  x  co-ordinate  of  the  point.  It  is 
obvious  that  for  the  diffuser  surface  the  variable  *  re¬ 
duces  itself  to  the  variable  x. 


It  is  evident  that  all  points  on  the  same  Mach  line 
(see  figure  10)  have  the  same  value  for  their  co-ordinates 
x  >  especially,  for  the  first  Mach  line,  it  is  x  =  0  ,  and 

therefore,  formula  (hi)  gives  immediately: 


51 


i 

u-(o/y  )  _  n  ~  i 

Uo  ~  / 


(U8) 


result  which  was  already  obtained,  for  a  more  general  case, 
in  section  5?  equation  (8). 

It  is  convenient,  by  analogy  with  the  function  0"(|)  , 
to  define  a  function: 


O'(f.v)  = 

u(0,y) 


(U9) 


where  \  ,  and  y  is  treated  as  a  parameter.  This 

function  is,  of  course,  a  generalization  of  the  function 
and  may  be  treated  as  a  limit  of  the  function  (com 
pare  equation  16): 


when  6  0  ,  and  y  =  const. 

From  (ii7)  and.  (Zj.8)  one  gets  immediately: 


5±S'1_ 


. .  l  r>T.  — * 


\  +  53+6y  1±9X  Tt_ 

5  2T6  > 


(50) 


As  one  should  expect,  in  the  limiting  case  y  =  l,the 
variable  reduces  itself  to  \  ,  and  we  have 
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O-(^l)  =  0{{) 

which  can  be  verified,  comparing  the  formulae  (50)  and  (i+O) 
For  y  =  oo  we  have  a  second  limiting  case.  For  small 
values  of  l  the  function  O'(f,oo)  is  determined  by  the 
series: 


=  l-fi  +&r- 


83 

204-9 


+ .... 


(5D 


and  for  arbitrary  \  : 


O'(l.oo) 


?  -It 

_  _i _ Mt)  e  dt 

VH’J  Kto+^tfCt)  t* 

0 


(51') 


The  shape  of  the  functions  0'(|,l)  and  (T(\,<x>) 
does  not  differ  essentially.  Both  functions  have  the  same 
initial  values  (=1)  and  both  approach  zero,  for  large  val¬ 
ues  of  I  ,  as  .  The  slope  of  O'  (|,d)  at  the  point 

i  =  0  is  -  ~  and  that  of  O'(jr(00)  is  -  -| 

The  functions  O'  (\,y)  for  y  =  const,  differ  still 
less  from  the  function  0'(],l)  =  CT(^)  ■  ,  because  their 
values  are  situated  between  these  two  limiting  cases  (i.e. 
y  =  1  and  y  =  co  ).  The  functions  cr ( ^ ,  d)  =  O’  (\)  ,  0(\,2.) 
and  O'(j.oo)  ,  for  values  of  I  in  the  range  between  0 
and  1,  are  shown  on  the  graph  9  • 
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19.  THE  ANGLE  OF  INCIDENCE 

In  the  case  of  axial  flow  without  an  angle  of  inci¬ 
dence.,  the  angle  between  the  flow  direction  and  the  dif¬ 
fuser  surface  is  constant  for  a  given  cross  section.  How¬ 
ever,  for  the  flow  with  an  angle  of  incidence  cc  this  is 
no  more  the  case.  For  this  reason  it  is  convenient  now  to 
assume  a  cylindrical  co-ordinate  system  x,r,  i?'  such  that 
(i)  -  the  x  axis  coi.cides  with  the  axis  of  the  dif¬ 
fuser 

(ii)  -  the  polar  angle  -fr  is  chosen  in  such  a  manner 
that  the  meridian  plane  which  contains  both 
the  direction  of  the  undisturbed  velocity,  and 
the  x  axis  corresponds  to  the  value  =  0  (see 
figure  11). 


Fie*  11 


5h 


It  is  obvious  that  the  angle  of  incidence  oc  of  the 
meridian  component  of  the  undisturbed  velocity  (for  a  merid 
ian  plane  with  an  angle  $  )  is: 

tahoc(^)  =  t.ana  cost?* 

(see  figure  12) 


or,  for  small  values  of  cc  : 

oC(t»  -  oc  Cos  $ 

Naturally,  the  component  of  the  undisturbed  velocity 
Uro  ,  normal  to  the  meridian  plane,  is  not  disturbed  by 
the  diffuser  and,  consequently,  does  not  influence  the  po¬ 
tential  'f  of  the  distrubed  portion  of  the  flow. 

As  is  known  from  the  theory  of  small  disturbances  for 
a  body  of  revolution,  the  part  of  the  potential  cp  which  is 
generated  by  the  disturbances,  in  the  case  of  a  flow  with 
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an  angle  of  incidence,  can  be.  decomposed  into  two  poten¬ 
tials; 

1)  the  first  one,  if*  ,  equal  to  the  corresponding  po 
tential  of  the  flow  without  an  angle  of  incidence, 
and 

2)  the  second  one,  or  cross  flow  potential, 

^  -  A  costf' 

where  A  is  to  be  obtained  from  the  boundary  condi¬ 
tion,  namely; 


yk  =  -  U<»«  cos  ft 

0  r 

Taking  into  consideration  that  the  operational  form 
of  is  already  known  (see  section  111)  and  can  be  repre¬ 
sented  as 

P  [<p*]  —  B  Ko  ((irp)  . 

(where  B  does  not  depend  on  r  ),  the  following  result  can 
be  immediately  obtained; 


PM=  P  [A  £?«*»] 

-  C  jip  C0ST>  Ki(prp) 
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where  £  is  the  value  to  be  obtained  from  the  boundary  con¬ 
dition  just  mentioned.  The  operational  form  of  this  condi¬ 
tion  is: 


.  3  r  J 


_  a 


£PM  - 


-  |-'Cp2  pZ  COST?'  [K0fprp)+'Kl(prp)]  - 


=  "  OtUoo  costf 


The  value  of  C,  resulting  from  the  last  formula,  is 


_ 2  oc  Uqq _ 

|3l  pJ[Ko((irp)  +  ^i(firp)] 


and  therefore: 


p  I  -  _  2alAoo  cost?'  Ki(prp)  (52) 

J  Pp[Kc(prp)  +  Ki  (jirp)] 


Recalling  that: 


and  therefore 


(x) 


All 

bx 
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P[u2(x)]  =  P[|?]  = 

one  obtains  finally: 

ulzM  _  T)"1  _  2<ifPr^  cos# 

oclXoo  _  P  [Ko(jir  p)  +  Kz(prf3')]_ 


(53) 


which  is  the  operational  representation  of  the  additional 
Uj.  velocity  generated  by  the  angle  of  incidence  of  the 
stream. 

The  non-operational  form  of  u.^-)  can  be  obtained 
by  a  procedure  analogous  to  that  of  section  l5»  It  may  be 
performed  by  changing,  in  a  suitable  way,  the  path  of  in¬ 
tegration  of  the  integral  involved  in  the  operation  P_1 
or  by  approximate  methods.,  such  as  finding  the  asymptotic 
expansions  of  the  Bessel  functions  Ko  ,  Ki  and  K  2.  for 
large  values  of  p 

It  is  also  possible  to  extend  these  results  to  the 
whole  space  around  the  diffuser,  and  for  a  non-limiting 
case  of  0  =  0,  by  use  of  the  same  methods  and  considera¬ 
tions  analogous  to  those  made  in  order  to  generalize  the 
results  of  the  case  without  an  incidence  angle,  as  shown 
in  the  former  sections. 

The  only  difference  in  comparison  with  the  former 


58 


methods  will  appear  because  of  the  singular  points  of  the 
functions  ^^(h)  and  K2f2‘)  .  In  view  of  these  singular 
points,  it  is  necessary  to  take  into  consideration  the  re¬ 
sidues  of  the  integrand,  while  changing  the  path  of  inte¬ 
gration  In  order  to  perform  the  operation  P'1  in  the 
formula  (53) • 
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Graph  1  -  u-velocity  along  a  straight  line  diffuser  ~  characteristics  method 


■raph  2  -  u-yeloe ity  along  a  straight  line  diff'usar  -  characteristics  method* 
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Graph  3  -  u-valocity  along  a  straight  line  diffuser  -  characteristics  method, 


Graph  5  -  Comparison  between  the  function  <T(|)  and  its  approximations 
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Graph  6  -  Comparison  between  the  characteristics-  method  and  the  analytical  solution 


Graph.  7  -  Comparison  between  the  characteristics  method  and  the  analytical  solution. 


Graph.  8  -  Comparison  between  the  oharacteristies  method  and  the  analytical  solution 


